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Extreme-mass-ratio inspirals (EMRIs) and intermediate-mass-ratio inspirals (IMRIs) — binaries in 
which a stellar-mass object spirals into a massive black hole or other massive, compact body — are 
important sources of gravitational waves for LISA and LIGO, respectively. Thorne has speculated 
that the waves from EMRIs and IMRIs encode, in principle, all the details of (i) the central body's 
spacetime geometry (metric), (ii) the tidal coupling (energy and angular momentum exchange) 
between the central body and orbiting object, and (iii) the evolving orbital elements. Fintan Ryan 
has given a first partial proof that this speculation is correct: Restricting himself to nearly circular, 
nearly equatorial orbits and ignoring tidal coupling, Ryan proved that the central body's metric is 
encoded in the waves. In this paper we generalize Ryan's theorem. Retaining Ryan's restriction 
to nearly circular and nearly equatorial orbits, and dropping the assumption of no tidal coupling, 
we prove that Thome's conjecture is nearly fully correct: the waves encode not only the central 
body's metric but also the evolving orbital elements and (in a sense slightly different from Thome's 
conjecture) the evolving tidal coupling. 

PACS numbers: 04.30.Db, 04.70.-s, 04.80.Nn 



I. INTRODUCTION AND SUMMARY 

The LIGO-GEO-VIRGO-TAMA network of broad- 
band ground-based laser interferometers, aimed at de- 
tecting gravitational waves in the high-frequency band 
10-10^ Hz, is already operating at or near its initial de- 
sign sensitivities. In the next decade, LISA (the Laser 
Interferometer Space Antenna) will open up the low- 
frequency gravitational- wave window (10~^-0.1 Hz). 

Among the most important sources of gravitational 
waves for LISA are extreme-mass-ratio inspirals (EM- 
RIs), which are systems in which a small object (with 
mass /i Mq) orbits a supermassive black hole or 
other central body (boson star or soliton star Q 

or naked singularity or ...) with mass M ~ lO^M©. Re- 
cently, Brown and collaborators Q have estimated that 
advanced detectors in LIGO (the Laser Interferomet- 
ric Gravitational- Wave Observatory) may detect up to 
~ 10 — 30yr~^ intermediate mass ratio inspirals (IMRIs), 
which are analogous to EMRIs but have less massive cen- 
tral bodies (masses M in the range of ^ 10^ — 

Thorne has conjectured "'^ that the waves from an EMRI 
or IMRI contain, encoded in themselves (at least in prin- 



^ Thome's conjecture has grown over time. Originally, in the early 
1990s, he conjectured (or, more precisely, asserted!) that the 
waves encode "a portion" of the spacetime geometry (e.g. p. 326 
of [^). By 1994, when Fintan Ryan proved his theorem, Thorne 
was arguing that the entire spacetime geometry would be en- 
coded (see, e.g., the introduction to Ryan's paper d). In 2002, 
when thinking about how LISA might test the laws of black-hole 
physics, Thorne realized that the tidal coupling might also be 
encoded along with the central body's spacetime geometry; see 
Ref. 0]. Only recently, when advising the authors about their 
research, did Thorne realize that the evolving orbital elements 
might also be extractable (private communication). 



ciple) : (i) the spacetime geometry (metric) of the massive 
central body, (ii) the tidal coupling (evolving rate of en- 
ergy and angular momentum exchange) between the or- 
biting object and the central body, and (iii) the evolving 
orbital elements. This conjecture (which has been par- 
tially proved; see below) has motivated placing EMRIs 
high on LISA'S list of target sources @, and has mo- 
tivated research to: (a) prove Thome's conjecture with 
the widest generality possible, or, if it is false, determine 
what information actually is encoded in the EMRI and 
IMRI waves (b) develop data analysis techni ques 

for searchingfor EMRI and IMRI waves in LISA [3 
and LIGO data; (c) scope out the accuracy with 
which LISA and LIGO can extract the encoded informa- 
tion from EMRI and IMRI waves (and if the central body 
appears to be a black hole, the accuracy with which its 
properties agree with those of a hole) [l3|| [l^: and (d) de- 
velop data analysis techniques for extracting the waves' 
information p^ . 

Fintan Ryan Q has proved a theorem that is an impor- 
tant step toward verifying Thome's conjecture. Specif- 
ically, he has proved that it is possible in principle to 
recover the full spacetime geometry from EMRI waves 
under the following assumptions: i) the central body is 
general-relativistic, stationary, axisymmetric, reflection- 
symmetric, and asymptotically-flat (SARSAF), ii) the 
small object travels on a nearly circular and nearly equa- 
torial orbit, and iii) there is no tidal coupling. Moreover, 
Ryan has shown that the multipole moments that deter- 
mine the spacetime geometry are redundantly encoded in 
the gravitational waves and can be extracted using ei- 
ther of the two precession frequencies (about a circular 
orbit and about the equatorial plane) or the waves' phase 
evolution. 

The purpose of this paper is to generalize Ryan's theo- 
rem. We retain assumptions (i) and (ii) (SARSAF space- 
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time and nearly circular, nearly equatorial orbit) but re- 
lax assumption (iii) by allowing for a small amount of 
tidal coupling. We show that in this case. Thome's con- 
jecture is nearly correct: the waves encode not only the 
central body's metric but also the evolving orbital ele- 
ments and (in a sense slightly different from Thome's 
conjecture) the evolving tidal coupling. (Assumption (ii), 
that the orbit is nearly circular and nearly equatorial, is 
relaxed in a companion paper by Li [^). 

Motivated by the result of Fang and Lovelace [1^1 that 
the only unambiguous part of the tidal coupling is the 
time-dependent, dissipative portion (at least when the 
central body is a non-spinning black hole and the orbit 
is large and circular), we characterize the tidal coupling 
by the rates of energy and angular momentum exchange 
between the central body and the orbiting object, i?body 
and ibody (Throughout this paper, a dot means deriva- 
tive with respect to the coordinate time t, which is the 
time measured by an inertial observer in the asymptot- 
ically flat region of the spacetime.) Actually, we only 
need to consider i?body, because once it is known, ibody 
can be deduced from the standard energy-angular mo- 
mentum relation for circular orbits and their influence 
on waves and tides, E — fiorbit^- (Here fiorbit is the or- 
bital angular velocity, which is the same as the waves' 
observed primary angular frequency aside from a factor 

. . . 

This paper is organized as follows: In Sec. |TT1 we begin 

by noting that, when there is a small amount of tidal 
coupling (as we assume), then the redundancy in Ryan's 
analysis is broken. One can still use Ryan's algorithm for 
the precession frequencies to recover the central body's 
spacetime geometry. Then, by making use of the ob- 
served (time-independent) spacetime geometry and the 
measured, evolving amplitudes associated with the pre- 
cession frequencies, one can also recover from the EMRI 
waves the evolving orbital parameters. Having relied on 
non-dissipative aspects of the waves to deduce the space- 
time geometry and orbit, one can then — as we show in 
Sec, inn — use the waves' dissipation-induced phase evolu- 
tion to deduce the tidal coupling. 

In our somewhat delicate discussion of deducing the 
tidal coupling fSec. lIIin . we begin by noting that the sum 
of the power radiated to infinity and the power fed into 
the central body via tidal coupling, i?totai = Eoo + i?body 
is equal to the power lost from the orbit, which can be 
deduced from the waves' observed phase evolution. The 
central body influences this observed -Etotai in two ways: 
(i) by generating a nonzero i?body, the quantity that in- 
terests us, and (ii) by very slightly altering Eao- To 
help quantify these two body influences, in Sec. lIIIBl we 
show how one can deduce, from the observations, the rate 
£'cx3NBi that energy would be radiated to infinity if there 
were no body influences. The difference between the mea- 
sured i?totai and the deduced ScxjNBI is the influence of 
the body's structure on the total energy loss from the or- 
bit, ii^totai,Bi = -Etotai " -EooNBi- This measured/dcduccd 
body influence on the total energy loss consists of two 



tiny pieces: the power that actually goes into the body 
via tidal coupling, ii^body, and the body's tiny influence 
on the power radiated to infinity, i?ooBi = ~ -BooNBi: 

-Etotal,BI = -BooBI + Ehody ■ (1) 

In principle (as described above), from the observa- 
tional data plus general-relativity theory we know the 
body's influence on the total energy loss ii^totai,Bi with 
complete precision. This is not quite what Thorne conjec- 
tured, but it is close, and it is the only complete-precision 
statement we have been able to make about measuring 
the influence of tidal coupling. 

Thorne conjectured we could deduce i?body from the 
observed waves. This, in fact appears not to be possible 
(in principle) with complete precision. However, we ar- 
gue in Sec. IIII CI and the Appendix that, if the central 
body is highly compact, then the unknown Eo^bi will be 
smaller than i^body by <C 1, where v is the orbital 

velocity and n is some high power; and we show that, 
when the body's external metric is that of Schwarzschild 
or Kerr, then n ~ 5. As a result, aside from a very small 
0(w") uncertainty due to the influence of the body on 
the energy radiated to infinity, the tidal coupling power 
^-body is equal to the known influence of the body on the 
total energy loss i^totai.Bi- 

A brief conclusion is made in Sec. |lVl 

II. EXTRACTING THE SPACETIME 
GEOMETRY AND ORBITAL ELEMENTS 

Aside from allowing tidal coupling, we treat the same 
class of EMRIs as did Ryan: 

First, we assume the central body's exterior spacetime 
is a vacuum, stationary, axisymmetric, reflection sym- 
metric, and and asymptotic flat (SARSAF) solution of 
Einstein's equations. The exterior spacetime metric can 
be written as (e.g., Eq. (7.1.22) of Ref. ^) 

ds^ =- -F{dt - ujd4)f + 4[e2'^(V + dz^) + p'^dcf)^], (2) 

where F,uj and 7 are functions of p and \z\. In SARSAF 
spacetimes, there is a one to one correspondence between 
the spacetime metric and a series of scalar multipole mo- 
ments {M2i,S2^+l), z = 0, 1, ■ • • [ii,[l^. Here Mq = M 
is the mass of the central body. Si is its spin, M2 is its 
mass quadrupole moment, etc. To extract the geome- 
try of the spacetime surrounding the central body, it is 
sufficient to extract the multipole moments {Me, Sg} 

Second, we let a small object with mass fi <^ M move 
about the central body in a nearly circular, nearly equa- 
torial orbit. 

For precisely circular, equatorial, geodesic motion, the 
waves obviously have a single fundamental frequency fJ^ 
that is associated with the circular motion (j) = fl^t. 
When the geodesic orbit is slightly nonradial, it is easy 
to show that its radius p undergoes periodic motion with 
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some angular frequency fip; and when slightly nonequa- 
torial, its vertical coordinate z undergoes periodic motion 
with another angular frequency fi^. These geodesic mo- 
tions give rise to gravitational waves that are triperiodic: 
a discrete spectrum with frequencies equal to Sl^, fJp, fJ^, 
and their sums and differences (see Ref. [sj for a proof, 
patterned after the proof by Drasco and Hughes [20] for 
the Kerr metric). The frequency difference fip — il^ shows 
up as an orbital periapsis precession, and il^ — fJ^ as an 
orbital plane precession; these precessions produce cor- 
responding modulations of the gravitational waveforms. 

In our case, the orbits are not geodesies; they evolve 
due to gravitational radiation reaction. Because of the 
extreme mass ratio, the radiation reaction can be de- 
scribed by the adiabatic approximation. In this approx- 
imation, on the timescalc of an orbital period, the small 
object moves very nearly along a geodesic of the cen- 
tral body's gravitational field. On a timescale much 
larger than the orbital period, the object moves from one 
geodesic to another as it loses energy and angular mo- 
mentum to gravitational radiation. It follows that the 
three frequencies {^^^(t), f2p(i), f2z(t)} each evolve with 
time on the radiation reaction timescale which is much 
longer than the orbital periods. 

In principle, a large amount of information can be en- 
coded in the time evolution of the waves' three funda- 
mental frequencies rip(t), ^zif'Yi and the complex 
amplitudes (amplitudes and phases) of the various spec- 
tral components. The largest amplitudes arc likely to 
be those for the second harmonic of fi^ and for the two 
precessions, /i2n^(i), hn^^n^{t), and hn^^n, (t)- We shaU 
call these the primary-frequency component, and the pre- 
cessional components of the waves. To simplify our prose, 
we shall refer to Op and as the "precession frequen- 
cies" even though the actual frequencies of precession are 
Op — O^ and Qz — O^. 

Thome's conjecture can be expressed mathematically 
as the claim that these time evolving frequencies and 
amplitudes encode fully and separably, 

1. the values of all the central body's multipolc mo- 
ments {Ml, Si}, 

2. the rates iSbody and ibody at which the orbiting 
object's tidal pull deposits energy and angular mo- 
mentum into the central body, and 

3. the time-evolving orbital elements, i.e. the orbit's 
scmi-latus rectum p{t), eccentricity e{t) and incli- 
nation angle t(t). 

Ryan's theorem @ states that, if there is no tidal cou- 
pling, then all the SARSAF moments {A'hi, 821+1} are 
encoded in the time evolving frequencies fully, separably, 
and redundantly. Ryan did not explicitly address the en- 
coding of the three orbital elements p{t), e{t) and t(i). 
However, their encoding is an almost trivial extension of 
his analysis: 

Specifically, Ryan noticed that the three fundamen- 
tal frequencies are independent of e and l to first order 



in these small quantities, i.e., they are functions solely of 
the moments and the semi-latus rectum p. One can elim- 
inate p by regarding the precession frequencies 0^ and Op 
as functions of the moments and O^, or equivalently as 
functions of the moments and the Post-Newtonian (PN) 
expansion parameter v = (A/O^)^/'^ 2± (orbital velocity). 
Expanding flz^v; Si, AIi) and Qp{v; Si, Mi) in powers of 
V, Ryan found the following pattern of coefficients (with 
each moment first appearing at a different power of v), 
from which all the moments can be extracted separably 
(Eqs. (18) - (19) of Ref. 0): 
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This result leads to Ryan's algorithm for extract- 
ing information. First, from the waves' observed time- 
evolving precession frequencies and time-evolving pri- 
mary frequency, one can deduce the functions O^_p(O0) 
and thence Oz^p(u); second, expanding in powers of v, 
one can then read out the multipole moments {Mi, Si} 
from either flz{v) or Op(u). 

We almost trivially augment onto Ryan's algorithm 
the following steps for extracting the time-evolving or- 
bital elements: Third, knowing the moments and thence 
the metric, one can use the geodesic equation to deduce 
p{t) from 0^(4). Fourth, one can use wave-generation 
theory and knowledge of the metric to deduce e{t) and 
L{t) from the amplitudes Hq^-ci^ and ft-o.-n^ of the wave 
modulations due to periapse precession and orbital plane 
precession. 



III. PROBING TIDAL COUPLING 

We now drop Ryan's restriction of no tidal coupling. 
This does not alter Eqs. ([3]) for Op and O^ as functions 
of V, i.e. of the orbital frequency O^, since all three fre- 
quencies only depend on the geodesic motion and hence 
only depend on the multipole moments {Mi, Si}. On the 
other hand, the evolution of the frequencies, as functions 
of time, will depend on the tidal coupling. 

More generally, we can divide the physical quantities 
of our analysis into two categories: i) "static:" those 
quantities related to the geodesic motion of the orbiting 
object, and ii) "dynamic:" those quantities related to 
the inspiral of the object (i.e., to the evolving rate at 
which the object moves from geodesic to geodesic). All 
static quantities are independent of tidal coupling and all 
dynamic quantities depend on it. 

This suggests that Ryan's analysis can be extended to 
include tidal coupling. First, the static quantities can be 
used to deduce the the central body's multipole moments, 
just as in Ryan's original argument as sketched above. 
Then, the dynamic quantities, combined with knowledge 
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of the spacetime metric, can be used to extract tidal- 
coupling information. This extension is discussed in the 
following subsections. 



A. The phase evolution when tidal coupling is 
neglected 

Following Ryan, we characterize the phase evolution of 
EMRI waves by the number of primary-frequency cycles 
of waves per logarithmic frequency interval, as a func- 
tion of the primary waves' slowly increasing frequency 
/ = ^^^/tt. This quantity can be written as (Eq. (4) of 
Ref. i) 



B. Tidal coupling and the phase evolution 

When tidal coupling effects are included, the redun- 
dancy is broken. The multipole moments {Mi, S^} can 
still [Eq. ([21)] be determined from ilp^z{^<j>), while (as the 
following discussion shows) , the tidal coupling can be de- 
termined from {Me,Se} and AN{f). 

As a preliminary to discussing this, we explain why it is 
sufficient, in analyzing tidal coupling, to focus on energy 
exchange between the orbit, the body and the waves, and 
ignore angular momentum exchange. Since the body is in 
a (nearly) circular, geodesic orbit, changes in its orbital 
energy and angular momentum are related by 
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by 



This phase evolution AN{f) can be measured 
gravitational-wave detectors with high precision. 

If there is no tidal coupling and no other influence 
of the structure of the central body on the waves, as 
Ryan assumed, then it is possible to read off the multipole 
moments (and also the small object's mass^) from a PN 
expansion of A7V(/) (Eq. (57) of Ref. 0): 
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Here "NBI" stands for no body influence and "TNILM" 
stands for terms nonlinear in lower moments. [Recall 
that V = {Mfl^y/^ = (7rA//)i/3.] So long as tidal 
coupling is negligible, then, the spacetime multipole mo- 
ments can be determined redundantly from either AN{f) 
[Eq. ([5])] or the periapse precession frequency Clp{n^) or 
the orbital-plane precession frequency ^^(fi^) [Eqs. ([3])]. 



^ The mass of the small object can be determined from AN{f) 
even when there is tidal coupling. The leading-PN-ordcr part 
of the energy flux (cquivalontly, the Icading-PN-ordcr part of 
AN(f)) is independent of tidal coupling. One can thus equate 
the leading-PN-order parts of A7V(/) and AATnbi [Eq. 10]. Af- 
ter inserting the mass M (obtained from one of the precession 
frequencies), one can solve for fi. The precession frequencies, in 
contrast, are independent of fi [Eqs. ((Sjl]. 



aside from second-order corrections due to the slight or- 
bital ellipticity and inclination angle. Our entire analysis 
is restricted to first-order corrections, so those second- 
order corrections are negligible. Similarly, since the en- 
ergy and angular momentum radiated to infinity are 
carried by the primary waves, with angular frequency 
Lu — 27r/ = 2il0 (aside from negligible contributions 
from the precessions, which are second order in the el- 
lipticity and inclination), each graviton carries an energy 
huj = 2Ml^ and an angular momentum 2h (with this last 
2 being the graviton spin). Therefore, the energy and 
angular momentum radiated to infinity are related by 



(6b) 



Conservation of energy and of angular momentum, to- 
gether with Eqs. ([5a|) and (j6bp . then imply that 
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(6c) 



for the energy and angular momentum deposited in the 
body by tidal coupling. Equations ^ imply that, once 
we understand, observationally, the energy exchange, an 
understanding of the angular momentum exchange will 
follow immediately. 

Now turn to the influence of the body's internal struc- 
ture on the observed energy exchange. 

The total rate that energy is lost from the orbit (which 
then goes to infinity and the body) is related to the phase 
evolution A^(/) by 



E, 



total 



■bit 



dEorhit df _ ^2 dEo-chit 1 



df dt 



df AN' 



(7) 



The phase evolution AA^ and the primary frequency / 
are known from observation, and, after using the preces- 
sion frequencies to compute the spacetime metric (Sec. 
HI)) , it is possible to compute dEd-hit/df via the geodesic 
equation'^. Thus everything on the right-hand side of 



To do this, first insert the multipole moments {Mi, Si} into 
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Eq. ([7]) can be determined from observed quantities, 
which means that i^totai is measurable. 

Another measurable quantity, we claim, is the rate that 
energy would be lost from the orbit if the body's structure 
had no influence. This quantity is [by analogy with Eq. 
0] 



total, NBI 



■bit 



1 



df AiVNBI 



(8) 



Knowing the moments as a function of frequency from 
measurements of the precessions, AiVNBi can be com- 
puted from the moments via Ryan's phasing relation'' ([5|). 
and, as we have seen, dEoihit /df can also be computed 
from the observations; so i?totai,NBi is, indeed, observ- 
able, as claimed. Therefore the influence of the body's 
structure on the orbit's total energy loss 
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total 



E, 



total, NBI 



(9) 



is also observable. 

This body influence on the total energy loss consists of 
two parts: the energy that goes into the body via tidal 
coupling, ii^body, and a tiny body-influenced modification 
of the rate that the waves carry energy to infinity 



£^total,BI — -Embody + -E^ooBI 



where 



E, 



ooBI 



— Eoo — £'total. 



NBI 



(10) 



(11) 



Thorne conjectured that the energy exchange due to 
tidal coupling, i^body, would be observable. We doubt 
very much that it is, since in general we see no way to 
determine the relative contributions of i?body and -EooBI 
to the observed total body influence i^totai.Bi- The best 
one can do, in general, in validating Thome's conjecture, 
is to extract the central body's total influence on the 
orbital energy loss, E'totai.Bi- However, in the special case 
of a body that is exceedingly compact, we can do better, 
as we shall explain in the next subsection. 



C. The dependence of Eoo on the central body's 
internal structure 

Consider a central body sufficiently compact that grav- 
ity near its surface blue shifts the orbiting object's tidal 



the geodesic equation. Then, solve the geodesic equation for 
the family of circular, equatorial orbits about the central body. 
Each orbit i will have a particular value of energy -Eorbit i a^nd 
frequency fi; this one-to-one mapping between Sorbit / can 
then be used to compute dEoi-bit/^f ■ 
* Ryan calculates the phasing relation to 2PN order (i.e., to 0(t!^) 
past leading order). By extending Ryan's calculation to higher 
post-Newtonian orders, the terms omitted from Eq. l(5j can be 
written explicitly. 



field, making it appear like ingoing gravitational waves 
as seen by stationary observers. This is the case, for ex- 
ample, when the central body is a black hole. Then, we 
claim, the ratio i?ooBi/£'body is very small: 



i^cxjBI 

^^body 



« 1 , 



(12) 



where n is a large number, very likely 5. For LISA, almost 
all of the wave cycles used in extracting information from 
the waves will be from radii where v < 0.5 so < 0.03. 
For example, for a Kerr black hole, if the spin parame- 
ter is a/M < 0.5, then at the innermost stable circular 
orbit, V < 0.5. Consequently, almost all of the measured 
^'totai.Bi will go into the body itself via tidal coupling, so 
^'body will be measured to good accuracy. 

To understand our claim that -EooBiZ-Ebody ~ v" for 
some large n, consider a central body whose external 
metric is that of a Kerr black hole. In this case, one 
can use the Teukolsky formalism [2]| (first-order pertur- 
bation theory in the mass ratio fi/M) to compute the 
energies radiated to infinity and tidally coupled into the 
central body. We have carried out that Teukolsky analy- 
sis for general a/M and present the details for the special 
case a = in the Appendix. Here we explain the under- 
lying physics. We begin with some preliminaries: 

We need only consider the primary-frequency waves, 
/ = fi^/TT, since they account for all the energy loss and 
transfer, up to corrections second order in the eccentricity 
e and inclination angle t. This means, correspondingly, 
that we can restrict ourselves to a precisely circular and 
equatorial orbit. The waves and tidal coupling then have 
predominantly spheroidal harmonic order £ = m = 2 
and frequency / (angular frequency lu = 27r/ ^ 2^1^). 
Since we only want to know, to within factors of order 
unity, the ratio i^coBi/^'body, it will be sufRcient to re- 
strict ourselves to these dominant £ = m ^ 2, uj = 2n^ 
perturbations. 

In the Teukolsky formalism, these perturbations are 
embodied in a radial "wave function" that can be nor- 
malized in a variety of different ways. The usual normal- 
ization, based on the Newman-Penrose field ip^, is bad 
for physical insight because it treats outgoing waves and 
ingoing waves quite differently; see the Appendix. One 
normalization that treats them on the same footing sets 
the radial wave function equal to that of the tidal gravi- 
tational field ( "electric-type" components of the Weyl or 
Riemann curvature tensor) measured by "zero-angular- 
momentum" observers, ZAMOs (a family of observers, 
each of whom resides at fixed radius r and polar angle 0) . 
We shall denote that tidal field [with e~"^* x (spheroidal 
harmonic) factored out so the field is complex, not real] 
by £. Another, closely related normalization for the ra- 
dial wave function sets its modulus squared equal to the 
rate of flow of energy. We shall denote this choice by \E'. 
At large radii, £ ~ {h+ + ihx) — w^(/i+ +«/ix ), where h+ 
and hx are the dimensionless gravitational wave flelds; so 
the radiated energy is Eoo ~ + |^ ~ {^/^)'^^?<n 
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FIG. 1: The renormalized tidal gravitational fields ^ pro- 
duced near a central body's surface and at large radii by the 
orbiting object, when the central body has the same exterior 
metric as a Kerr black hole. 



which tcUs us that "^oo ~ {'r/^)£oo- Near the body's 
surface (i.e. near where the horizon would be if the 
body were a Kerr black hole), the energy flux is _E ~ 
(r/a;)^|a^f P, where a is the Kerr-metric lapse function, 
which goes to zero at the horizon radius. (The ZAMOs' 
divergently large outward speed, relative to infalling ob- 
servers, causes them to see a divergently large tidal field; 
the factor corrects for that divergence; see, e.g., the 
discussion in Sec. VI. C. 2 of [l^].) Thus, in order to en- 
sure that the power flow is the square of the renormalized 
radial wave function, 



E - l^-l 



(13) 



we must renormalize the ZAMO-measured tidal field £ 
by 



~ {r/uj)£ at r ^ oo , 
4' ^ {o?r lijj)£ near body 



(14) 



With these preliminaries finished, we can give our 
physical argument for Eq. p2p in terms of the radial wave 
function 5*. Our argument relies on Fig. [T] 

If the central body is a Kerr black hole, then the 
boundary condition on 4* at its surface (the horizon) is 
purely downgoing waves, and at infinity, purely outgoing 
waves. The ratio of downgoing power at the horizon to 
out goin g power at infinity has the standard Kerr values 
[23t]23|: i^body/^^ooNBi ~ if the hole's spin angular 
velocity f2// is much less than the orbital angular veloc- 
ity Jl^; and £'body/£'ooNBi ^ if fin > 51^. (Here we 
have used the no-body-influence notation EqcNBI for the 
outgoing power because a central black hole's internal 
structure is unable to influence the waves radiated to in- 
finity.) Correspondingly, by virtue of Eq. (|13|) . the ratio 
of the downgoing field at the horizon to the outgoing 
field at infinity ^'ooNBI is 



V 
7,5/2 



for 



(15) 



This suppression of the downgoing field relative to the 
outgoing is due, mathematically, to a reflective effective 
potential in the wave equation that ^ satisfies (Fig. [T]). 
Physically, it is due to coupling of the field 5" to the 
central body's spacetime curvature. 



Now suppose the central body is not a black hole, but 
some other object so compact that its surface is well be- 
neath the peak of the effective potential. This mathemat- 
ical assumption is equivalent to our physical assumption 
that the ZAMOs see the downgoing field so strongly 
blue shifted by the central body's gravity that it looks 
like radiation. The only way, then, that the central body 
can influence the energy radiated to infinity is to reflect a 
portion of this downgoing radiation back upward. Math- 
ematically, this corresponds to replacing the black hole's 
downgoing boundary condition by 



(16) 



at some chosen radius just above the body's surface. Here 
and are the downgoing and upgoing components 
of ^f; see Fig. [TJ For simplicity, we shall assume that 
the amplitude reflection coefficient TZ is small, \Ti\ <^ 1. 
Otherwise we would have to deal with a possible resonant 
buildup of energy between the reflective central body and 
the reflective effective potential — though that would 
not change our final answer (see, e.g., the more detailed 
analysis in the Appendix). 

The upgoing waves have great difficulty getting 
through the effective potential. The fraction of the up- 
going power that gets transmitted through, successfully, 
is ^ (Muj)^ if the hole rotates slowly, and ^ {Mui)^ if 
rapidly [Eq. (8.83) of [1^ with £ = 2 and a^o = w]. Since 
the fields ^ are the square roots of the powers (aside 
from complex phase) and since Moj = 2M^^ = 2v^ , this 
power transmissivity corresponds to 



V 

,15/2 



for 



(17) 



Combining Eqs. (fT7|) . (fTBl) . and (fT51) . we see that 



^ r 



;BI 



)NBI 



for 



Oh > 00 



(18) 



If these two complex outgoing fields are not pre- 
cisely out of phase with each other (phase difference 
±7r/2), then the outgoing power is IvPooNBI + 'J'ooBiP — 
^ooNBiP + 2K (\l/ooNBi^^Bi)' which means that the ra- 
tio of the radiated body-influenced power to radiated no- 
body-influence power is 



)BI 



E^ 



)NBI 



,13 
,10 



for 



Oh < O^ 
Oh > Orf 



(19) 



In the unlikely case (which we shall ignore) that the two 
fields are precisely out of phase, the ratio will be the 
square of this. 

By combining Eq. with the square of Eq. (|15p . 

we obtain the ratio of the body-influence power radiated 
to infinity over the tidal coupling power into the central 
body: 



-EooBlZ-Ebody 



(20) 



7 



independent of whether the body rotates slowly or 
rapidly. This is the claimed result. 

If the central body's external metric is not Kerr, then 
the first-order perturbation equations for the orbiting 
body's spacetime curvature will probably not be sepa- 
rable in {r,9}, so the analysis will be much more com- 
plex. Nevertheless the physical situation presumably is 
unchanged in this sense: The body's spacetime curvature 
will couple to the perturbation field in such a way as to 
resist energy flow through the region between the body's 
surface and the object's orbit. Correspondingly, the per- 
turbation fields and power flows are very likely to behave 
in the same manner as for the Kerr metric, with the same 
final result, i?body/-EooBi ~ v" with n very likely still 5 
but possibly some other number significantly larger than 
one. 

If this is, indeed, the case, then for any sufficiently 
compact central body the power tidally deposited into 
the body Sbody will be very nearly equal to Etoted. Bii 
which is measurable; and therefore the tidal power will 
be measurable. 



review the standard Teukolsky perturbation formalism. 
Details can be found in e.g. Ref. [2^. To shorten our 
expressions, in this appendix we restrict ourselves to a 
nonrotating central body with external metric the same 
as a Schwarzschild black hole but with a finite reflectivity. 
The generalization to the Kerr metric is straightforward 
but with more cumbersome algebra. We have carried it 
out, obtaining the same result as is found by the physical 
argument in the text. 

In the Teukolsky formalism, people usually calculate 
the perturbation to a Newman-Penrose quantity that 
is related to the ZAMO-measured tidal field by a linear 
transformation of the basis vectors. This ipi can be de- 
composed into Fourier-Harmonic components according 
to 



V'4--/ diuY,R-^in^{r)_2Yi^{9,^)e-'^\ (Al) 

Irn 



IV. CONCLUSION 

In this paper, we have extended Ryan's analysis to 
show that in principle it is possible to recover not only 
the spacetime geometry of the central body, but also the 
evolving orbital parameters of the inspiraling object and 
the evolving tidal coupling between the small object and 
the central body. Therefore, in principle we can obtain a 
full description of the SARSAF spacetime, the tidal cou- 
pling, and the inspiral orbit from EMRI or IMRI wave- 
forms. In practice, the method of extracting the infor- 
mation is likely to be quite different from the algorithm 
we have presented here. 

Further generalizations of Ryan's theorem and devel- 
opment of practical methods to implement it are topics 
of our ongoing research. 
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APPENDIX A: AN EXPLICIT DERIV ATION OF 
RESULTS IN SEC, unci 



where -2Yim{0, 4>) the spin- weighted spherical har- 

il function i?, 
mogeneous Teukolsky equation 



monies. The radial function R,oim{r) satisfies the inho- 



r^a^—r~2{r-M)— + U{r 
dr 



dr 



Ru,im{r)^T^iU^2) 



where = I — 2M/r is the lapse function for the 
Schwarzschild metric. The expressions for the potential 
U (r) and the source T^im can be found in e. g. , Ref. [11], 
Eqs. (2.3), (Al). 

In order to solve this equation, we construct two lin- 
early independent solutions to the homogeneous Teukol- 
sky equation, which satisfy the following boundary con- 
ditions. 



4 ^ 4 iojr* 



IN 
wlm 



I \4 4 

(wr) a e 
{ujrY 



2M 



ulrn 



3 /^out ^iujr* 



}lm 



UP 



{LOT) a r^i^e + r^i^e , r 



2M 



3 iior* 



(A3) 



1. Teukolsky Perturbation Formalism 

In this subsection, we use the Teukolsky perturbation 
theory to justify our results in Sec. IIII CI We first briefly 



where d/dr* = a^d/dr. From these two homogeneous 
solutions, we can construct the inhomogeneous solution 
according to 



8 



Wronskian[i?X,fi™J 



(r) 



2M 



dr'R]^Ur')r.Ur')\,m 



r 



where 7^;m(r) = Tt^;,„(r)(r^ — 2Mr)~^. This solution has 
only outgoing waves at infinity and satisfies the purely 
ingoing boundary condition: (Ref. [23| . Eqs. (2.8) and 
(2.9)) 



Rulm{r ^ oo 
i?^/™(r ^ 2Af) 



\ 2 ^7 



2 rrlN Jiiujr* 



dm, 



re 



3 ^UP „4„4„-iwr* 



r a e 



(A5) 



where 



,IN,UP 



(A6) 



X I dr 

I2M 



{r'-2Mr)-'RZ-^^{r)T^,„,{r) 



At infinity, where the spacetime is almost flat, ?/'4 is di- 
rectly related to the outgoing gravitational wave strains 
according to 



V'4 



ihs 



(A7) 



and we can obtain the luminosity formula (Ref. [23| . 
Eq. (2.21)) 



(A8) 



2. Inner Boundary Condition 

The above purely ingoing boundary condition makes 
sense when the central body is a black hole because we 
know everything is absorbed at the horizon of the black 
hole. If the central body is some other kind of object, 
the only way it can influence the perturbation field Rujim 
just above its surface is by producing an outgoing-wave 
component via some effective reflectivity TZ. The result 
will be a modified field 



Ru^imir -> 2M) ~ e- 



(something)e*' 



(A9) 



The "something" will be proportional to TZ, and it will 
also have a peculiar radial dependence because V'4 relics 
for its definition on an ingoing null tetrad and thereby 
treats ingoing and outgoing waves in very different man- 
ners. 



functions were introduced and used in Sec. IIII CI the 
ZAMO-measured tidal field £ and a field 5" whose mod- 
ulus squared is the power flow, for both outgoing and 
ingoing waves. Those choices are good for Sec lIIICf s 
physical, order-of-magnitude arguments, but at general 
radii r they not related in any simple way to '04- A 
choice that is simply related to ip4 is the Regge- Wheeler 
function X, and we shall use it here. 

The radial wave function R for the Newman-Pcnrose 
-04 is related to the Regge- Wheeler function X by the 
Chandrasekhar transform, Eq. (A6) of Ref. [2^. This 
Chandrasekhar transform takes the form 



R: 



IN, UP 
ilm 



IN,UP^ yIN,UP 



where 



16(1 - 2iMuj){l ~ AiMuj){l + AiMuj) 



{l-l)l{l + l){l + 2)-12iMuj 



(AlO) 
(All) 



IN, UP 
An 



are 



is a second order differential operator, and X 
two linearly-independent solutions of the homogeneous 
Regge- Wheeler equation 



where 



dr* 



V{r) 



V{r) 



X^hnir) = 0, 



(A12) 



1(1-1) 



6M 



IN, UP 



The asymptotic expressions for X^^^ 
Eq. (2.7)) 



(A13) 
are (Ref. [H, 



IN 



X, 



ujlm 




2M 



A'. 



out — iojr 
ujlrn 



nil! pi' 



(A14) 



2M 



-oo 



Here we note that by the conservation of the Wronskian, 
it is straightforward to show that S'"-""' = A'n.out^ 



4. Eoo with a Reflective Inner Boundary Condition 



3. Chandrasekhar Transform 

To learn what the "something" should be, we can 
transform to a new radial wave function that treats ingo- 
ing and outgoing waves on the same footing. Two such 



Because the Regge- Wheeler function treats outgoing 
and ingoing waves on the same footing, the desired, re- 
flective inner boundary condition for it takes the simple 
form 



XX(r-2M) 



(A15) 
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Here is a new homogeneous solution of the Regge- 

Wheeler equation. 

This new homogeneous solution is a superposition of 
both ingoing and outgoing waves at the horizon. It is 
shown in Ref. [2^ that because the Regge- Wheeler func- 
tion treats outgoing and ingoing waves in the same man- 
ner, |7^p has the physical meaning of the energy flux 
reflectivity, i.e., the ratio between outgoing and ingoing 
energy flux at the horizon. 

The homogeneous solution l|A15[) which satisfies the 
new inner boundary condition can be constructed from 
the old homogeneous solutions: 



IN 



(A16) 



where 



From these new Z^^^ the calculation of the luminosity 
at infinity is straightforward. 

In Ref. [2^ Poisson and Sasaki have already worked 
out all the relevant formulae, so we only give the results. 
For the original expressions in Ref. [23| . please refer to 
Eq. (3.25) for A"\ Eq. (A7) for x™, x"^, Eqs. (5.4), 
(5.6), (5.11), (5.12) for 

The leading luminosity correction comes from the I = 
2, m = ±2 mode, and we have 



E rso = En 



Schwarzschild 



1 - 



128^^7^^;l3 



15/3i 



(A20) 



/3i = 1 



■RA 



out 

iolm 



/32 = 



n 



A" 



(A17) 



After doing an inverse Chandrasekhar transform, we 
obtain the corresponding homogeneous solution of the 
homogeneous Tcukolsky equation 



i?; 



IN 



IN 
!l7n 



R. 



UP 



(A18) 



Now we can replace i?™ by ^™ in Eq. (jXI)) to ob- 
tain the solution Rujim{T) which satisfies the inhomoge- 
neous Teukolsky equation with upgoing and downgoing 
waves at the horizon and purely outgoing waves at infin- 
ity. From this Ruhnir) we identify the new amplitudes 
Z'^im as in Eq. jM]): 



where v is the same FN expansion parameter as that in 
Sec. nil CI Unless the reflection coefficient TZ is precisely 
real, this gives 



E „a — Ea 



Schwarzschild 



256 , 



,13 



(A21) 



in agreement with Eq. (|19[1 . The change in -©body should 
be 



Ehody — -^body 



Schwarzschild \^ l/^ll"^ ^ 



^IN 



Z, 



IN 



P2X 

01 X, 



IN 

UP ^Lolm 
Lolm 



(A19) where /3i is defined in Eq. (|A17|) . 
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